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CLASSIFICATION  OF  FLOWS  IN  AIR  SHOCK  TUBES 


G.  M.  Arutyunyan 
(Moscow) 


Examined  are  the  possible  types  of  flow  with  constant  pressure 
behind  the  front  of  a  shock  wave  with  arbitrary  variation  of  initial 
pressures  in  the  sections  of  the  shock  tube  and  of  the  geometric  di¬ 
mensions  of  the  tube.  It  Is  assumed  that  the  low  pressure  section  of 
the  shock  tube  has  an  open  end. 

1.  Formulation  of  Problem.  Consider  a  shock  tube  of  constant 
cross  section  (Fig.  I)  having  a  low  pressure  section  with  an  open  end. 
Let  the  air  In  the  high  pressure  section  of  length  l  be  in  a  state  of 
compression  (c);  and  that  In  the  low  pressure  section  of  length  z,  In 
the  state  (0).  We  shall  assume  that  the  air  in  both  sections  Is  In 
thermal  equilibrium  with  the  ambient  medium. 

When  the  diaphragm  is  swiftly  demolished  (Fig.  2),  a  shock  wave 
(S)  begins  to  propagate  to  the  right  (through  the  low  pressure  medium) 
and  a  centered  rarefaction  wave  (i)  to  the  left  (through  the  high 
pressure  medium).  In  addition,  there  is  formed  contact  discontinuity 
K  which  Is  carried  to  the  right  by  the  flow.  At  this  discontinuity 
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all  the  parameters  of  the  medium  with  the  exception  of  the  pressure 
and  particle  velocity  undergo  jump  variations. 

The  pressure  p  in  the  front  formed  by  the  shock  wave  may  be 
determined  from  the  relationship  (cp.  Grib  [l]) 


yr2k[lLk~i)ptlpa  +  (k  —  l)\ 


2* 


i 


(1.1) 


where  pfi  and  p0  are  the  initial  pressures  in  the  sections  of  the  shock 
tube,  and  k  is  the  isentropic  exponent. 

After  the  raref action  wave  (i)  reaches  the  closed  end  of  the 
tube,  its  reflection  begins,  forming  a  new  wave,  the  front  of  which  is 
propagated,  to  the  right  with  the  velocity  of  weak  disturbances  u  +  c 
(u  is  the  particle  velocity,  and  c_  is  the  local  speed  of  sound). 

At  the  moment  the  shock  wave  reaches  the  open  end  of  the  shock 
tube,  a  rarefaction  wave  is  generated,  the  front  of  which  will  propa¬ 
gate  to  the  left  with  a  velocity  u  -  c.  • 


Fig.  3. 


Thus  in  the  general  case  two  disturbances  will  be  propagated  in 
a  shock  tube  with  respect  to  the  homogeneous  states  (s)  and  (K)  in 
opposite  directions:  a  disturbance  from  the  closed  end  of  the  shock 
tube  and  a  disturbance  from  its  open  end. 

It  is  known  that  at  a  fixed  distance  from  the  diaphragm  x 
(x  >  0)  the  variation  of  the  pressure  in  the  shock  tube  through  time 
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has  the  form  'shown  in  Fig.  3-.  The  shock  wave  has  some  constant  pres¬ 
sure  area  of  length  At,  which  arises  when  the  shock  wave  front  reaches 
the  given  point  and  disappears  when  the  disturbance  from  the  open  or 
closed  end  of  the  shock  tube  reaches  this  point. 

In  practical  work  on  shock  tubes  it  is  important. to  know  the  size 
of  this  area,  within  the  confines  of  which  the  homogeneous  flows  (S)  - 
and.  (K) ■ and. the  contact  discontinuity  dividing  them  are  located.  .  , 

In  general,  the  quantity  At  is  a  function  of  five  parameters: 
ps,  p0,  L,  V,  and  x. 

■  2.  ,  '.  'Some  Necessary  Relationships.  Let  us'  determine  the  conditions 

under  which  (Fig.  4)  the  right  hand  boundary  of  rarefaction  wave  (l) 
will:  propagate  to  the  left,  remain  in  place,  or  propagate  to  the 
right.  In  other  words,  it  is  necessary  to  find  the  conditions  under 
which 


“.<«*•  ;  “«> c* 


(2.1) 


It  may  be-  shown  that 


(2.2) 


Hence  the  conditions  in  (2.1)  may  be  presented  respectively  in 
the  form: 


Pt>(  2  ytt-  T  2  v 

pr  \k-\-i)  9  pc .  -}-  iy  .  »  pc 


u 

Jr  1 


(2-3) 


Converting  from,  p  to  p0  with  the  aid  of  (l.l),  we  get  for  the 
case  of  air  (k  =1.4) 


-^.<2.69.  . 

Pa 


Pa 


-  2.89, 


—  >2.89 


(2.4) 


respectively. 
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It  may  be  shown  [2]  that 


i  3-*  irt  - 


(2.5) 


After  solving  the  equations  of  thq  characteristic  7  and  the 
trajectory  of  the  contact  discontinuity  simultaneously  (Fig.  4)  for 
the  coordinate  of  the  point  £  we  obtain 

’  .  /  l2'6) 


From  the  simultaneous  solution  of  the  equations  of  the  character¬ 
istic  71  and,  the  trajectory  of  .the  front  of  the  shock  wave,  the  fol¬ 
lowing  expression  may  be  obtained  for  the  x-eoordinate  of  their  point 
of  intersection  d 


2/P  c‘  j  \  1 
1“.-®  ek\*k) 


1  3"* 


(2.7) 


where  D  is  the  velocity  of  the  front  of  the  shock  wave.  Let  us  intro¬ 
duce  the  dimensionless  quantity 


20 


u.-i-e—D  c_ 


1  rt 


(2.8) 


Then  it  can  be  maintained  that  the  disturbance  from  the  closed  end  of 
the  shock  tube  is  successful  in  leaving  the  front  of  the  shock  wave 

behind  over  the  length  L,  if  L/ 1  >  a>x;  is  not  successful  in  leaving 

the  front  of  the  shock  wave  behind  over  the  length  L,  if  L/l  <  cdx. 

After  solving  the  equation  of  the  characteristic  72  (Fig.  5) 

together  with  the  equation  of  the  trajectory  of  the  contact  discon¬ 
tinuity,  we  obtain  the  following  expression  for  the  coordinate  of  the 
point  fz. 


Ws-M  ,  _HD  +  e.-uJ 

cJD  .  I  ej) 


(2.9) 


Introducing  next  the  dimensionless  -quantity 


i  :rk 


2 D  4v  /  \  ;  k~i 

“1=r  U^rct  —  ut  ck  {  ek  ) 


(2.10) 


-and  making  use  of  (2.6)  and  (2.9),  we  come  to  the  following  conclu-  .  , 
slon:  the  disturbance  from  the  closed  end  of  the  shock  tube  encounters 
the  surface  of  the  contact  discontinuity  before  the  disturbance  from 
the  open  end,  if  L /l  >  o>2;  the  .  disturbance  from  the  open  end  of.  the 
shock  tube  encounters  the  surface  of  the  contact  discontinuity  before 
the  disturbance  from  the  closed  end,  if  L/l  <  o>2;  both  disturbances 
encounter  the  surface  of.  the  contact : discontinuity  simultaneously,  if 
L/l  =  cd2. 

Simultaneous  solution  of  the  equations  of  the  characteristics  7 
and  73  allows  us' to  obtain  for  the.  x-coordinate  of  the  point  h  . the 
expression  ■'  ,  * 


..  **  Vvj  . 


(2.11) 


Then;  Introducing  the  dimensionless  quantity 


(Oj  =. 


2C'D(ct-U') 

k*  +  “,)  —  «,)  (  «*./ 


1  ark 


(2.12) 


It- is  possible  to  assert  that 


X/|<°  when  L> 1  < *h>0  nheni/l> 


<*1 


The  intersection  point  e_- of  the  .  characteristics  71  and  73  has  the 
coordinate 


UC—L!D)  +  L  (ut  +  ct)  —  xc  (u, —  e,) 
__ _  -  ; 


(2.13) 


Let'  us  introduce  the  dimensionless  quantities 

(c,\  Tp, 

ck  [2 c,D  —  (u.  +  ck)  {D+C'—  «,)]  (  ck  ) 

,  nr..  .  \  1  *"*  1  rk 

c.n(Ul-ck)  fcUlFi  _  u,  —  ek  I  c„  WiFi 


Oi  = 


a,  = 


/c»\TFI  _  l  Cp  '\7 

V* (P  +  V  j  ’  3  ct  (  et  j 


(2.14) 


Jo- 


Then 


Z  when  Z//< at,  xh  =  Z  when  L/l—a xh<[Z  wi*n  Z//>a, 

■  xb'>xi  when  L  / 1  <C. nj, .  x6  =  xf  when  Z //  =  a2,  •  x6<Xy  when  Z//>o» 

xt>Z  wfrei  Z/2<a3,  xb  =  Z  when  L/l  —  a 3.  x6<^Z  when  £/l>n. 


■For  classification  of  the  flows  in  the  case  u  >  c_  it  is  advan- 

S  5 

tageous  to  introduce  the  dimensionless  quantities 

l  r»~*  ’ 


■  •  ,Pi=«o.;  =  f-r)2*'1-.  3»=“* 

.  e*  \ck  J 


(2.15) 


Then  , 


xtt'^L  when  Z /.f  <CPi> 
x£>Z  when  Z  /  /<p2. 
.xt>  Z.  '  when  Z  /  2  <  Ph. ' 


xd  =  L  when  Z//  =  P>.  ■  ;  xd <  Z  .  when  £/*!> Pi 

xc=  Z  itoi  ,,Z/ 2  =  p2,  ■  .  xc<;Z  When  Z./1>P» 

xb— L  to  Z/I=  Pi,.  tt<i  itei  l;l>Pi 


Obviously,  in  order  to  determine  the  size  of  the  constant  pres-- 
sure  area '  it  Is  necessary  to  single,  out  the  region  of'  constant  pres¬ 
sure  p  =  pg,  formed  by  the  regions'  of  the  homogeneous  flows  (S')  and  • 
(K)'  in  the  xt  plane..  /'  ..  •  '  '' 

",  •  In  general  the  boundaries  of  the  constant  pressure  region  are 
made  up  of  the  trajectory . of  the  front  of  the  shock  wave'  and  the  char¬ 
acteristics  7P;'7,.  7i,-  72>  and  j3,  and  for  this  reason  a  constant 
pressure  region  of  a  particular  configuration  in '  the  xt  plane  jrill 
hereafter  be  called  a  characteristic  regime. 

Let  us  show  how  many  different  regimes  there  can  be  in  a  shock 
tube  for  arbitrary  variation -of . the  parameters  p  ,  p<>,  L,  and  l  and 
what  their  distinguishing  characteristics  are. 

5.'  Possible  Characteristic  Regimes  for  the  Case  ug  <.  e^.  It  is 
possible  . to  show  that  for  subsonic:  flows  the  inequality 


><i)j 


(5.1) 


will  always  hold  behind  the  contact  discontinuity. 
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The  case  ug  <  c^.  We  shall  choose  the  regime  depicted  In  Fig.  5- 
Xt  is  characterized  "by  the  failure  of  the  disturbance  from  the  closed 
end  of  the  tube  ,  to  leave  the  front  of  the  shock  wave  behind  over  the 
distance  L;  the  disturbance  from  the  open  end  of  the  shock  tube  en¬ 
counters  the  surface  of  the  contact  discontinuity  before  the  disturb¬ 
ance  from  its  closed  end;  x^- <  0.  In  other  words,  this  regime,  if  it 
■is  possible,  should  be  characterized  by  simultaneous  fulfillment  of 
the  relationships:  L/t  <  L/l  <  ui2,  h/l  <  u>3.  Taking  into  account 
inequalities  (3. l),  "we,,  shall  require  that  L/l  satisfy  the  condition 

.  •  .  '  .  '•  If  I  J  (5*2) 

Obviously,  here’ the  first  two  conditions  are  fulfilled  automat¬ 
ically.  Consequently,  this,  characteristic  regime  is  possible  in  prin¬ 
ciple,  for  which  it-  is  enough  to  require  fulfillment  of  condition. 

:  (3*2).  .t  '-  • 

With  the  aid  of  the  equations  of  the  trajectory  of  the  shock 
wave  .front  and  the  equations  of  the  characteristics  72  and  73  it  is 
pqssible  to  determine  the  size  of  the  constant  pressure  area  in  the 
whole  interval.  0  <  x  <  L2..  Other  possible  types  of  characteristic 
regimes  for  the  case  u  <  c^.  are  presented  in  Fig.  6.  Also  stated 
there  are  the  conditions  under  which  they  occur. 

The  case  u?  =  c^.  The  possible  characteristic  regimes  for  the 
case  of  sonic  flows  behind  the  contact  discontinuity  are  presented  in 
Fig-  7* 

4.  ,  Possible  Characteristic  Regimes  for  the  Case  .c^  <  ug  <  cg. 

We'  note  that  when  the  flow  behind  the  contact  discontinuity  in  a  shock 
tube  becomes  supersonic,  it  still  remains  subsonic  behind  the  shock 
wave  front.  The  possible  chai’acteristic  regimes  for  the  case 


c.  <  u  <  c  are  considered  below  taking  into  account  the  fact  that 
iC  s  s 

under  these  conditions  the  relationship 

Mr->  <*>,  >  Oj  >  a»  ■'  (4.1) 

Is  always  true. 

As  an  example,  let  us  consider  the  regime  depicted  .in  Fig.  8.  . 

It  is  characterized  by  the  fact  that  the  disturbance  from  the  closed 
end  of  the  shock  tube  does  not. have  time  to  overtake  the  shock  wave 
.front  on  length  L;  the  disturbance  from  the  open  end  of  the  shock  tube 
encounters  the  surface  of  the  contact  discontinuity  before  the  dis¬ 
turbance  from  the  closed  end: 

*h>-L>  *b>xr  '  *b>L. 

.  .  In  other  words,  this  regime,  if  it  is ' possible,  should  be  char¬ 
acterized  by  simultaneous  fulfillment  of  the. relationships: 

A /:<<■>,.  ///<«-.  L/ I<o,.  £/./<  ns.  i/t<a, 

¥e  shall  require  that  Jj/l  satisfy  the  condition  h/l  <  03. 
Then  the  other  conditions  will  be  fulfilled  automatically  on  the  basis 
of  inequalities  (4.1).  Consequently,'  this  characteristic  regime  Is 
possible  in  principle. 

Other  characteristic  regimes  possible  for  the  case  c^  <  ug  <  eg 
are  presented  in  Fig.  9  with. a  statement  of  the  conditions  under  which 
they  occur. 

Thus,  when  c.  <  u  <  c_  there  can  be  no  more  than  II  character- 
istic  regimes  In  the  shock  tube;  the  actual  type  of  flow  is  determined 
by  the  relationship  between,  the  dimensionless  quantities  cdx,  cd2,  ax, 
02,  and  03,  and  the  value  of  h/l. 
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It  Is 


5-  Possible  Characteristic  Regimes  for  the  Case  u  >  c  . 
_ _ _ s  —  s 

possible  to  show  that  the  relationship 

Pi  >?=>?..  •  (5.1) 

is  always  true  for  supersonic  flows  behind  the  shock  wave  front. 

Cqnsider  the  regime  depicted  in  Fig.  10.  It  is  characterized  by 
the  relationships:  xd  >  L,  xc  >  L,.  and  x^  >  L.  In  other  words,  in  -  / 
order  to  realize  such  a  flow  the  conditions  L/l  <  px,  L /l  <  p2,  and 
L/l  <  p3  must  be  fulfilled  simultaneously. 

‘  We  shall  require  that  the  inequality  L/l  <  p3  be  fulfilled.  Then 
the  others  will  be  fulfilled  automatically  on  the  basis  of  (5-l).  Con¬ 
sequently,  this  characteristic  regime  Is  also  possible  In  principle. 
Other  possible  types  of  flow  when  u  >  c  ■  are  given  In  Fig.  II. 

Thus,  when  ug  ;>  eg  the  characteristic  regime  Is  determined  by 
the  relationship  between  L/l  and  the  dimensionless  quantities  px,  p2, 
and  p3. 

6.  Conclusion.  It  has  been  established  above  that  for  arbitrary 

variation  of  the  parameters  p  ,  p0,  L,  and  l  there  may  exist  in  an 

s 

air  shock  tube  12  characteristic  regimes  when  u  <  c.  ,  11  when 

S  K 

ck  <  us  <  es^  and  7  when  ug  >  cg.  Altogether,  therefore,  no  more 
than  thirty  types  of  flow  may  exist  in  the  shock  tube. 

It  can  be  shown  that  to  assume  the  existence  of  any  other  char¬ 
acteristic  regimes  In  addition  to  the  thirty  presented  above  would 
mean  violating  conditions  (5.1),  (4.1),  and  (5.1). 

The  dimensionless  quantities  0)3.,  ca2,  0)3,  ax,  a2,  a3,  px,  p2,  and 
p3  are  functions  only  of  the  ratio  pg/p0.  Presented  In  Fig.  12  are 
the  curves  of  these  functions,  which  divide  the  graph  plane  Into  a 
number  of  regions  corresponding  to  the  thirty  possible  characteristic 
regimes  considered  above.  This  representation  Is  convenient  for 
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determining  tine  type  of  flow  for  given  values  of  the  parameters  ,p  , 

s 

Po,  L,  and  l  (or  pc,  p0,  L,  and  z).  In  fact,  for  this  purpose  it  is 
enough  to  form  the  ratios  p  /p0 .  :and  L/z  and  then  establish  in  which 
of  the  characteristic  regions  in  Fig.  ±2.  the  point  determined  by  these 
ratios  in  the  graph  plane  belongs. 


The  author  expresses  his  deep  gratitude  to  Ya.  B.  Zel’dovich, 

A.  S.  Kompaneyets,  and  to  Kh.  A.  Rakhmatulin  for  discussing  this  work 
and  for  their  valuable  advice." 
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